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ABSTRACT 


A computer  program  called  TRIBRID  is  described  which  can  rapidly 
calculate  the  smooth  water  resistances  of  some  advanced  bibrid  and 
tribrid  ship  concepts  that  are  supported  by  buoyant,  dynamic,  and 
powered  aerostatic  forces.  These  forces  are  produced  by  one  or  more 
hull-strut  combinations  similar  to  a demi  SWATH  in  conjuction  with 
large  hydrofoils  and/or  air  cushions.  For  a given  displacement  and 
shape,  the  size  of  each  concept  is  determined  and  various  drag  components 
along  with  some  interference  drag  are  then  calculated  for  a calm  water, 
fixed  zero  trim  and  free  to  heave  operating  condition. 


INTRODUCTION 


The  primary  intent  of  this  report  is  to  describe  the  TRIBRID 
computer  program  which  can  quickly  calculate  the  smooth  water 
resistance  of  some  of  the  advanced  hybrid  ship  concepts 
investigated  by  J.  Meyer  (Reference  1)  and  D.  Jewell  and  of  many 
possible  and  foreseeable  modifications  of  these  concepts.  The  terms 
hybrid  and  tribrid  refer  to  a ship  which  is  supported  by  a combination 
of  buoyant,  hydrofoil,  and  powered  air  cushion  generated  lift.  Bibrid 
vehicles  have  two  of  the  above  support  modes  and  tribrids  have  three. 

The  computer  program  was  initially  developed  for  the  Hydrofoil 
Small  Waterplane  Area  Ship  (HYSWAS)  and  has  been  expanded  to  handle 
the  Small  Waterplane  Area  Air  Cushion  Ship  (SWAACS)  and  many 
variations  of  the  Small  Waterplane  Area  Twin  Hull  (SWATH)  ship  with 
air  cushion  and/or  hydrofoil  support.  In  addition,  the  program  may  be 
used,  with  some  modification,  for  predicting  some  of  the  calm  water  drag 
components  of  the  Hydrofoil  Air  Cushion  Ship  (HYACS)  and  the  Large  Hybrid 
Hydrofoil  Ship  (LAHHS).  Figure  1 shows  some  of  these  concepts. 
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TRIBRID  PROGRAM  DESCRIPTION 


The  TRIBRID  program  was  written  to  quickly  estimate  the  resistance 
of  hybrid  and  tribrid  ship  concepts  which  depend  on  buoyant,  dynamic 
or  air  cushion  sustention  devices.  The  buoyant  displacement  is  derived 

from  a slender  body  of  revolution  (called  a hull  1n  this  report)  operating 
under  the  free  surface  and  also  from  a low  aspect  ratio  strut  situated  on 
top  of  the  body  and  piercing  the  free  surface.  The  dynamic  lift  is  due 
to  very  large  subcavitating  foils  moving  through  the  water,  and  the  air 
cushion  lift  is  derived  from  pressure  supplied  by  pov/ered  fans  aboard 
the  ship.  The  hull-strut  combination  can  be  used  singly,  in  pairs  as  in 
a catamaran  or  in  triplets  as  in  a trimaran.  The  foils  may  be  cantilevered 
from  the  hull  or  strut  like  the  wings  of  an  airplane  or  they  may  span 
two  adjacent  hulls  or  struts  depending  on  the  desired  arrangement.  The 
air  cushions  utilize  the  struts  as  side  seals  and  they  may  be  placed 
between  adjacent  struts.  The  resistance  of  a large  mumber  of  possible 
hull-strut,  foil  and  cushion  arrangements  can  be  evaluated  by  the  TRIBRID 
program.  A few  of  these  arrangements  are  shown  in  the  TRIBRID  program 
input  section  of  this  report. 

Existing  equations  and  computer  programs  which  were  available  at  the 
inception  of  the  project  were  gathered  to  predict  the  drag  components 
associated  with  the  foils,  air  cushions  and  hull-strut  combinations  and  were 
incorporated  into  the  TRIBRID  program.  In  the  area  of  foil  wave  drag  and 
the  wave  drag  of  flared  strut*:  no  suitable  existing  program  was  found. 


A subroutine  was  written  to  evaluate  the  wave  drag  of  large  foils 
operating  at  shallow  submergences  and  an  existing  computer  program 
was  modified  to  take  into  account  the  effect  of  strut  flare  on  the 
wavemaking  drag  of  the  hull -strut  combination.  Since  it  was  not  the 
intent  of  this  project  to  develop  new  technology,  a number  of 
foreseeable  interference  effects  between  the  foils,  hull -strut  combination, 
and  air  cushions  along  with  some  other  drag  components  have  not  been 
evaluated.  A complete  list  of  anticipated  drag  components  that  are 
not  evaluated  is  found  in  the  section  on  the  limitations  of  the  TRIBRID 
program. 

Since  some  drag  components  are  not  accounted  for  by  the  TRIBRID 
computer  program,  it  is  inadvisable  to  make  resistance  comparisons 
between  a concept  evaluated  by  the  TRIBRID  program  and  an  already 
existing  ship  whose  resistance  is  known.  The  TRIBRID  program  neglects 
certain  resistances  and  therefore  the  comparison  of  the  resistance 
generated  by  TRIBRID  to  a known  ship  resistance  must  be  considered 
speculative  to  a certain  extent. 

Approximate  methods  are  used  (See  the  Appendices  for  details) 
for  the  computation  of  the  various  drag  components.  Hard  data  are  not 
available  to  verify  the  validity  of  most  of  these  formulations.  It  is 
thus  conceivable  that  the  ranking  of  widely  different  concepts  based 
solely  on  these  powering-performance  computations  may  not  be  born  out 
in  reality.  For  example,  if  one  compares  a concept  having  mainly 
dynamic-lift  sustension  with  another  that  is  mainly  cushion-supported. 
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an  error  in  the  relative  ranking  could  be  caused  by  the  overpredicted 
cushion  drag.  Additional  error  in  the  ranking  can  be  introduced  if 
non-optimum  versions  of  the  concepts  are  used  to  represent  each  component. 

It  should  be  emphasized  here  that  the  TRIBRID  program  deals 
only  with  resistance  and  powering.  Powering  is  addressed  superficially, 
and  only  because  it  is  necessary  to  calculate  the  total  required 
power  for  valid  comparisons  between  cushion  and  non  cushion  concepts. 

An  assumed  fan  efficiency,  duct  efficiency,  propulsive  coefficient  and 
the  calculated  effective  power  and  fan  ideal  power  are  used  to  determine 
the  total  required  power. 

The  feasibility  of  a configuration  is  not  examined  with  regard 
to  propulsor  sizing  or  propulsor  efficiency  calculations,  nor  are  any 
basic  Naval  Architectural  feasibility  questions  with  regard  to  static 
stability,  motions,  manuevering,  structures,  machinery  arrangement, 
fuel  and  payloads  investigated. 
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TRIBRID  PROGRAM  INPUT 


The  required  input  to  run  the  TRIBRID  program  consists  of  the 
following  data  cards.  Numbers  in  [brackets]  are  SI  unit  conversion 
factors,  (i.e.,  quantity  in  British  Unit  x Conversion  Factor  = 
quantity  expressed  in  SI  units.) 


Card  1 12A6  Format 

Columns  1-72  Enter  concept  name  or  title 

Card  2 8F10.6  Format  - Constants 

Columns  1-10  RHO  - density  of  water  slugs/ft^  [515.3788  K^/m^] 

2 2 

11-20  GRAV  - gravitational  constant  ft/sec  [.3048  m/s  ] 

21-30  PI  - Value  of  pi  3.14159... 

31-40  E - Logarithmic  base  2.71828... 

41-50  GAMMA  - Euler's  constant  .55721... 

51-60  Vise.  - Kinematic  Viscosity  of  water  ft/sec^  [10.76423  m/s^] 

61-70  RHOAIR  - Air  density  sluqs  [515.3788  K„/m^] 

T 9 

ft'^ 

Card  3 8F10.6  Format 

Columns  1-10  DELTCF  - Correlation  coefficient  (.0005  for  example) 

11-20  FKS  - Equivalent  sand  grain  diameter-used  as  a measure  of 
foil  roughness  - inches  [2.54  cm] 

Card  4 8F10.6  Format 

Columns  1-10  DELTA  - Total  buoyant  displacement  at  rest.  Tons  Weight 

[RHO  X GRAV  *0.9910  m^] 
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11-20 

21-30 

31-^0 

41-50 

51-60 


61-70 

Card  5 8F10.6 

Columns  1-10 

11-20 

21-30 

31-40 

41-50 

51-60 


CWP  - Waterplane  coefficient  of  the  strut 

CP  - Hull  prismatic  coefficient 

Cl  - Strut  inertial  coefficient  (see  appendix  B for 
definition) 

EXP  - Strut  flare  exponent  enter  1.0,  2.0,  or  3.0  only. 

LCBOLH  - Ratio  of  the  hull  longitudinal  center  of 
buoyancy  measured  from  midships,  LCB,  to 
the  hull  lencth,  LH  (Positive=LCB  fwd  of  ^ ) 

LCFOLS  - Ratio  of  the  strut  longitudinal  center  of  floatation 
measured  from  amidships,  LCF,  to  the  strut  length,  LS 
(Positive=LCF  forward  of  (JJ  ) 

Format  (see  Figure  2 for  a schematic  showing  of  the 
geometric  quantities  discussed  below) 

TZODH  - Ratio  of  the  strut  thickness,  TZ,  at  strut  bottom 
to  horizontal  hull  diameter,  DH 

TR  - Ratio  of  the  strut  thickness,  BSH,  at  the  top  to  the 
thickness,  TZ,  at  bottom 

HSODV  - Ratio  of  strut  height,  HS,  to  the  vertical  hull 
diameter,  DV 

DVODH  - Ratio  of  hull  vertical  dimater,  DV,  to  horizontal 
diameter,  DH 

LSOLH  - Ratio  of  strut  length,  LS,  to  hull  length,  LH 

DHOLH  - Ratio  of  horizontal  hull  diameter,  DH,  to  hull  length,  LH 


Card  6 81 

Column  5 
in 
15 
20 
25 
30 


Format  - Concept  Configuration  Card 

NS  - Number  of  struts 

NSH  - Number  of  hull  strut  intersections 

NH  - Number  of  hulls 

NCUSS  - Number  of  cushions 

NF  - Number  of  foil  sets.  (NF=2  for  tandem  foils) 
NFC  - Number  of  foil  sections 
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A foil  set  may  consist  of  one,  two,  or  three  sections.  A one  section 
foil  set  is  a foil  spanning  two  hulls.  A two  section  foil  set 
consists  of  two  outboard  wing  foils  on  a single  hull,  or  on  two  hulls. 
A three  section  foil  consists  of  a center  spanning  foil  betv/een  two 
hulls  and  two  outboard  wing  foils.  See  Figure  3 for  some  of  the 
various  possible  foil,  hull,  strut  and  cushion  combinations. 

Card  7 8F10.6  Format  - VKNOT(I) 

Enter  the  eight  desired  speeds  in  knots.  [.5144444  m/s]  They  may  be 
the  same  speed  or  different  speeds. 


Card  8 8F10.6  Format  - CPSF(I)  [47.88026  Pa] 

2 

Enter  the  eight  cushion  pressures  in  pounds/ft  corresponding  to  the  above  speeds. 
Enter  zeroes  if  the  concept  does  not  use  aerostatic  lift. 


Card  9 8F10.6 

Columns  1-10 


Format 

LCOLS  - Enter  the  ratio  of  cushion  length  to 
strut  length 


(1.0) 


11-20  SLOLC  - Enter  the  ratio  of  the  combined  forward  (0.5) 
and  aft  seal  length  to  cushion  length 

21-30  CD  - Enter  aerodynamic  drag  coefficient  for  the  (0.3) 
craft 

31-40  HEOLS  - Enter  the  ratio  of  the  average  super-  (.06) 
structure  height,  HE,  to  strut  length 


Card  10  8F10.6  Format  - SODH(I)  ft.  [.3048  m]  - See  Figure  2 


Enter  eight  ratios  of  the  centerline  separation  distance  SEP,  between 
adjacent  hulls  to  the  horizontal  hull  diameter,  DH.  For  single 
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hull -strut  concepts  enter  zeroes.  For  catamaran  and  trimaran  type 
concepts  SODH  must  be  greater  than  1.0.  For  trimaran  concepts,  SODH 
is  the  ratio  of  the  distance  between  the  centerlines  of  the  center 
hull  and  the  side  hull  to  the  horizontal  hull  diameter. 


Card  11  8F10.6  Format  - CSTRUT(I)  ft.  [.3048  m] 


Enter  eight  values  of  the  longitudinal  distance  in  ft.  between  the 
centerline  of  strut  and  centerline  of  hull.  (+  is  fwd  of  hull  center- 
line) 


Card  12  8F10.6  Format 


Columns  1-10  FAC  - Foil  average  chord  ft.  [.3048  m] 

11-20  FATC  - Foil  average  thickness  to  average  chord  ratio 

21-30  V/TB  - Foil  wing  tip  to  body  distance  ft.  [.3048  m] 

if  NFC  on  card  6 is  1 then  set  WTB  equal  to  zero. 


31-40  AZ  - Foil  lift  curve  slope.  Change  in  foil  lift 
coefficient  per  angle  of  attack  in  degrees 

41-50  FOIST  - Tranverse  distance  in  ft  [,3048  m]  between 
the  tip  of  a central  foil  section  and  the 
body  centerline.  For  a central  foil  at  the  bow 
or  stern,  FOIST  is  zero. 

51-60  FVDIST  - Vertical  distance  in  ft  [.3048  ml  between 
the  centerplane  of  a central  foil  and/or 
wing  foil  and  the  hull  axis. 


Card  13  Format  810.6  -FLOAD(I) 

Enter  eight  values  of  foil  loading  corresponding  to  the  eight  speeds  on 
card  7.  Units  are  lb  [47.88026  Pa]  Enter  zeroes  if  the  craft  has  no 
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Card  14 


8F10.6  Format  - ETAPC 


Enter  8 propulsive  coefficients  corresponding  to  each  of  the  speeds  on 
card-  7.  ETAPC  must  be  greater  than  0. 

Card  15  8F10.6  Format  - ETAF 

Enter  eight  fan  efficiencies  corresponding  to  each  of  the  speeds  on 
card  7.  ETAF  must  be  greater  than  0.0. 

Card  16  8F10.6  Format  - ETAD 

Enter  eight  duct  efficiencies  corresponding  to  each  of  the  speeds  on 
card  7.  ETAD  must  be  greater  than  0.0. 

Card  17  815  Format 

Columns  1-15  NSTART  - Enter  99999  if  another  concept  is  being  run 

and  insert  cards  1-17  for  new  concept.  Enter 
zero  for  the  last  case  or  if  only  one  case  is 
run. 
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SUMMARY  OF  CALCULATED  LIFT  & DRAG  COMPONENTS 

The  following  is  a list  of  the  lift  and  drag  components  evaluated 
by  the  TRIBRID  program. 

HULL  - consists  of  a slender  body  of  revolution  below  the  free  surface. 
The  calculated  lift  and  drag  components  are; 

A)  Buoyant  lift 

B)  Friction  drag 

C)  Eddy  drag 

D)  Wave  making  drag 

E)  Wave  making  interference  drag  between  hull(s)  and  struts(s) 

STRUT  - The  strut  connects  the  upper  superstructure  or  cross-structure 
of  a concept  with  the  hull.  It  is  long  and  slender  and  may  be 
wall  sided  or  flared. 

A)  Buoyant  lift 

B)  Friction  drag 

C)  Eddy  and/or  tip  drag  (Tip  drag  is  calculated  only  for 

concepts  that  have  no  hull  under  one 
or  more  of  the  struts) 

D)  Wave  making  drag 

E)  Wave  making  interference  drag  between  strut(s)  and  hull(s) 

FOIL  - The  foil  has  a subcavitating  section  shape  and  is  attached  to 
the  hull  or  strut  like  a wing.  The  foil  can  also  span  two 


hulls  or  struts. 


A)  Dynamic  lift 

B)  Buoyant  lift 

C)  Friction  drag  (See  Appendix  A for  the  calculation  of  the 

foil  friction  coefficient  which  is  expected 
to  be  high  due  to  roughness  and  fouling 
resulting  from  long  periods  of  foil  submergence) 

D)  Induced  drag 

E)  Wave  making  drag 

F)  Eddy  making  drag 

AIR  CUSHION  - The  air  cushion  is  contained  between  two  adjacent  struts 
and  fore  and  aft  seals. 

A)  Cushion  aerostatic  lift 

B)  Wave  making  drag 

C)  Seal  drag 

D)  Momentum  drag 

E)  Fan  ideal  power  (calculated  for  purposes  of  comparing  cushion 

and  non  cushion  craft) 

SUPERSTRUCTURE  - The  aerodynamic  pressure  drag  is  calculated. 

The  equations  for  calculating  the  above  listed  lift  and  drag 
components  are  presented  in  appendices  A,  B,  and  C.  These  appendices 
correspond  to  the  modular  subroutines  in  the  TRIBRID  computer  program 
which  evaluate  foil  system  performance,  hull -strut  performance  and 
cushion  system  performance.  Each  of  the  appendices  has  its  own  list 
of  symbols  identical  to  the  symbols  in  the  original  reference  sources. 
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The  subroutine  corresponding  to  each  appendix  is  designed  so  that 
it  can  easily  be  replaced  by  a new  subroutine  describing  the  per- 
formance of  a different  lift  device  that  might  be  employed  on  a 
future  concept.  In  fact,  the  present  TRIBRID  nrogram  was  developed 
this  way  from  a January  1975  computer  program  that  had  the  capability 
of  evaluating  concepts  only  with  buoyant  and  dynamic  support  and  no 
air  cushion  support. 


TRIBRID  PROGRAM  OUTPUT 


Three  pages  of  output  are  generated  by  the  TRIBRID  program. 

The  first  contains  Quantities  which  do  not  vary  with  speed  such  as 
displacement,  hull,  strut,  foil  and  cushion  dimensions,  foil  lift 
curve  slope  and  number  of  hulls,  struts,  foils  and  cushions.  The 
various  hull,  strut,  foil,  cushion  and  superstructure  drag  components 
for  maximum  of  eight  specified  input  speeds  are  printed  in  tabular  form 
on  the  second  page  in  English  Units  and  on  the  third  page  in  SI  Units. 
The  tables  also  include  the  following. 


a)  Assumed  propulsive  coefficient,  fan  efficiency,  and  duct 
efficiency.  These  efficiencies  are  needed  in  order  to  compare 
the  total  power  required  by  craft  with  air  cushion  support  to 
craft  without  air  cushion. 

b)  Buoyant,  dynamic  and  aerostatic  lift  expressed  as  a percentage 
of  the  total  craft  weight.  The  percentages  are 'calculated  in 
order  to  place  each  concept  on  the  sustension  triangle  proposed 
by  Jewel  (Reference  2), 

c)  Cushion  pressure  - specified  for  each  speed. 

d)  Strut  wetted  surface,  strut  draft  on  the  cushion  side,  and  the 
total  draft  - these  are  all  quantities  which  change  with  the 
heave  of  the  craft. 

e)  Fan  flow  rate  - depends  on  cushion  pressure  iS  well  as  on  speed. 

f)  Foil  loading  - pressure  prescribed  by  the  user. 

g)  Foil  angle  of  attack  and  foil  lift  coefficient.  The  program  can 
generate  foil  angles  beyond  stall  if  the  foil  loading  specified 
is  too  high  for  a particular  speed.  The  user  should  check  that 
the  foil  angle  and  foil  lift  coefficient  are  below  the  stall 
values.  If  not,  the  loading  should  be  decreased. 

h)  Sum  of  the  drag  components  calculated 

i)  Effective  power 
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j)  Fan  ideal  power 

k)  Total  shaft  power  - the  total  shaft  power  is  the  sum  of  the 
effective  power  divided  by  the  propulsive  coefficient  plus  the  fan 
ideal  power  divided  by  the  product  of  the  fan  and  duct  efficiences, 

A sample  computer  output  is  shown  in  Figure  5.  The  sample  is 
for  a twin-hull, foil  and  cushion  supported  concept,  the  midship 
section  of  which  is  shown  in  Figure  4.  The  sample  output  is  included 
for  illustrative  purposes  only  and  not  for  making  any  evaluation  of 
the  concept. 


FIG.  4 MIDSHIP  SECTION  SCHEMATIC  CORRESPONDING  TO  SAMPLE  OUTPUT 
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FIGURE  5 - Computer  Output  for  Sample  Case 


TWIN  HULL  TRinRID  WITH  FOIL  AND  CUSHION  SUPPORT  txAMPLE  CASE 

velocity  5.14  7.73  10.29  12.86  15.43  18.01  20.58  23.15  M/S 

ASSUMED  PROPULSIVE  COEFFICIENT  .<;00  .500  .500  .500  .500  .500  .500  .500 

ASSUMED  FAN  EFFICIENCY  .500  .500  .500  .500  .500  .500  .500  .500 

ASSUMED  DUCT  EFFICIENCY  .700  .700  .700  .700  .700  .700  .700  .700 
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FIGURE  5 - Computer  Output  for  Sample  Case 


LIMITATIONS  OF  THE  TRIBRID  PROGRAM 


The  following  lift  and  drag  components  are  not  evaluated  by  the 
TRIBRID  program: 

1)  Any  wavemaking  interference  drag  between 

a)  foil  and  hull 

b)  foil  and  strut 

c)  foil  and  cushion 

d)  cushion  and  hull 

e)  cushion  and  strut 

2)  Any  boundary  layer  type  interference  drag,  for  example;  any  change 
in  foil  lift  or  drag  due  to  operating  a part  of  the  foil  in  the  hull 
boundary  layer  is  not  accounted  for. 

3)  Drag  due  to  any  part  of  the  ship  above  the  struts  coming  in  contact 
with  the  water.  This  is  especially  important  in  the  low  speed  drag 
prediction  for  the  HYSWAS  concept  which  relies  on  the  upper  hull 
for  static  roll  stability  prior  to  take-off.  The  upper  hull  must 
be  in  the  water  until  sufficient  speed  is  attained  for  the  foils  to 
provide  roll  stability.  The  resistance  of  the  upper  hull  is  not 
evaluated;  the  magnitude  of  expected  HYSWAS  hump  drag  at  take-off  will 
depend  upon  the  design  of  the  upper  hull  and  the  proportions  of 
buoyancy  and  dynamic  lift  selected  for  a particular  design. 

4)  Any  drag  or  lift  due  to  cavitation  and/or  ventilation  of  the  foils, 
strut  or  any  other  part  of  the  ship.  Fully  wetted  flow  is  assumed. 

5)  Drag  due  to  a non  zero  trim,  yaw  or  roll  angle.  The  TRIBRID  program 
assumes  that  the  craft  has  zero  trim,  zero  yaw  and  zero  heel  attitude. 
No  moment  or  equilibrium  calculations  are  performed. 
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6)  Drag  due  to  the  sinkage  of  the  hull  caused  by  dynamic  forces 

7)  Lift  and  Drag  generated  by  control  surfaces  such  as  rudders 
or  small  foils., 

8)  Appendage  drag  due  to  shafting,  bosses,  shaft  supports  or  other 
appendages 

9)  Rough  water  drag.  Only  smooth  water  performance  is  considered. 

10)  Spray  drag 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  TRIBRID  computer  program  is  a useful  tool  for  quickly  calculating 
the  resistances  of  many  advanced  ship  concepts.  The  program  is  limited 
in  that  no  hard  full  scale  data  exist  to  verify  some  of  the  drag 
component  predictions.  In  addition,  certain  forseeable  interference 
drags  are  neglected  because  there  is  no  accurate  way  of  calculating 
them.  It  is  strongly  recommended  that  resistance  experiments  be 
performed  as  a concept  in  order  to  predict  its  resistance  with  greater 
confidence. 

Additional  research  and  development  is  needed  to  improve  the 
powering  prediction  capability  of  the  TRIBRID  program.  At  present 
there  is  no  acceptably  accurate  tool  for  estimating  the  drag  of  the 
following  drag  components:  hull-foil-cushion  wavemaking  and  viscous 

interference  drag,  the  drag  due  to  cavitation  and/or  ventilation  of 
the  foils  and  struts,  spray  drag  of  the  strut  and  the  wavemaking  drag 
of  the  unconventional  upper  hulls. 
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APPENDIX  A 

FOIL  SYSTEM  PERFORMANCE 


k 


LIST  OF  SYMBOLS  FOR  APPENDIX  A 


AR 

B(k) 


D 

Dr 


"4 

E(k) 

On(^) 

K(k) 


R 

R 


cav 


S 

U 


2 

aspect  ratio  of  hydrofoil  = (2b)  /s 
modified  complete  elliptic  integral 
B(k)=l/2[E(k)-(l-k)^K(k)] 

foil  friction  coefficient  based  on  total  wetted  surface 

total  foil  drag  =D.+D.+D^+Dc 

foil  eddy  making  drag 

foil  friction  drag 

foil  induced  wave 

foil  wave  making  drag 

complete  elliptic  integral  of  the  second  kind  of  argument  (k) 
Bessel  function  of  the  first  kind  of  order  n 

complete  elliptic  integral  of  the  first  kind  of  argument  (k) 

total  foil  lift  =L  +AL,+aL^+L„ 

0 1 2 B 

lift  due  to  static  buoyant  displacement  of  foil 

aerodynamic  value  of  foil  lift 

g/U^ 

reynolds  number  based  on  average  foil  chord 

plan  area  of  the  foil 
free  stream  velocity 
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b 

“^avg 

c(y) 

g 

h 

t 

avg 

t(y) 

6 

r(y) 

c 

p 

0 

A 


Y 

AL 

AL, 


absolute  angle  of  attack,  measured  from  the  free  stream 
to  the  zero-lift  direction 

slope  of  the  lift  coefficient  curve 
half  span  of  the  foil 

average  chord  length  of  elliptic  planform  foil 
chord  length  of  the  foil  c(o)=c^ 

gravitational  acceleration  constant 

depth  of  submergence  of  foil  measured  from  trailing  edge 
to  mean  free  surface 

average  thickness  of  elliptic  planform  foil 

thickness  of  foil  t(o)=t^ 

Froude  number  with  respect  to  b 
distribution  of  the  circulation  strength  r(o)=r^ 
elevatio'-’,  of  the  disturbed  free  surface 
water  density 


Froude  number  with  respect  to  h= 


JL 

gh 


1 


depth  - 1/2  span  ratio  =h/b=B/a 
Euler's  constant  0.5772... 
lift  do  to  presence  of  free  surface 
lift  due  to  wave  effects 


APPENDIX  A - FOIL  SYSTEM  PERFORMANCE 


The  foil  prediction  method  is  based  upon  lifting  line  theory  with 
linearized  free  surface  assumptions  as  outlined  in  Reference  3 which 
considers  a hydrofoil  of  finite  span  moving  with  constant  velocity  at  a 
fixed  distance  below  the  surface.  See  Figure  lA. 


Fig.  lA-  Coordinate  system  for  the  hydrofoil  motion. 


The  details  of  the  general  formulation  of  the  problem,  the  boundary  con- 
ditions used  and  the  solution  of  the  Laplace  equation  are  presented  in 
Reference  3 and  will  not  be  repeated  here. 


All  the  following  equations  for  foil  performance  are  applicable  to 


a foil  with  elliptical  distribution  of  the  circulation  strength  r as 
given  by  equation  lA.  The  equation  numbers  in  brackets  are  the  equation 
numbers  as  they  appear  in  Reference  3. 


r(y)  = r 1 - ^ , |yUb, 
° 1 b'^ 


r(y)  = 0 , 


|y|  2 b. 


The  lift  and  drag  equations  are  presented  in  terms  of  a depth  Froude  number 
o,  a span  Froude  number  e and  a 1/2  span-depth  ratio  A,  as  shown  in 
equation  2A  and  3A 


The  equations  for  lift  are: 


where 


L = + aL^  + aL^  + Lg 


pC  ( 


- 3tiA  \ 


3^5  + a' 


1 - 3a"  jlog  2 - I + J log 


+ h(x  loqA)  as  A-»-0 


[78a] 


a4  = !!o!„l,  f,.^ 

4^  i ^ 


1 (t  + I09  f)  (1  + 0(X^  1)  [80]  (7A) 


r 


is  the  aerodynamic  value  for  a foil  in  an  infinite  medium.  L-j  is 
independent  of  gravitational  effects  and  represents  the  lift  due  to  the 
mean  free  surface.  I2  is  the  lift  caused  by  surface  wave  effects.  Lg 
represents  the  static  buoyant  foil  lift  and  is  calculated  from  the 
displaced  foil  volume.  The  foil  is  assumed  to  be  elliptic  in  planform 
with  a constant  thickness  to  chord  ratio  as  shown  in  Figure  2A 


The  drag  equations  are: 


D=D^  +D2+D2+D^+D^+Dp 


TrpF 


” 2 
J,  (pb) 


f,  - 4 u 


TTpr 


r CO 


D2+D3  = 4 


irpr 


(llA) 

CM 

C 

0. 

[60] 

(12A) 

J,^(ub) 

[61] 

(13A) 

- xJi+i^ 

IT  V 


1 


/i  • -2 


+ X-/ 


^ / 


V2  2 _ 2 


— sec  6 
a 


J,(lseA  sine)  ■;  ^^^5,  (UA) 

2 ! F63T 

sec  0 sine  ^ ■' 


is  the  induced  drag  due  to  the  trailing  vortices  of  the  hydrofoil. 

The  combination  of  D2’*'^3  '"spesents  the  total  contribution  of  the  mean 
free  surface  effect  which  favorably  decreases  the  total  drag  especially 
when  h is  small.  As  h tends  to  zero,  D2^'^3  cancel  Is  . D4  represents 
the  wave  drag  due  to  the  downstream  wave  formation  which  results  from  the 
gravity  effect.  In  equation  13A,K(k)  and  E(k)  denote  the  complete  elliptic 
integral  of  the  first  and  second  kind  respectively  and  are  evaluated  by 
the  NSRDC  Mathematics  Laboratory  library  function  CELLI.  The  integral 
portion  of  equation  14A  is  evaluated  using  standard  numerical  techniques. 


Dp  represents  the  viscous  friction  drag  and  is  given  by  equation  15A. 


“f  ' ^ ‘avg  '>) 

^avg  is  ^he  average  chord  length  and  the  quantity  (4  c^^^  b)  , spresents  the 
total  wetted  surface  of  the  foils. 

is  the  foil  friction  drag  coefficient  for  hydraulically  smooth 
foils.C^  is  given  by  the  Schoenherr  friction  line  as  given  in  equation 
16A  Reference  6 


0.242 


(16A) 


where  is  the  Reynolds  number  based  on  the  average  foil  chord  length. 

Coi  V 

For  foils  which  are  not  hydraulically  smooth,  the  foil  friction  coefficient 
is  given  by  equation  17A  Reference  6 


C 


f 


(1.89  + 1.62  log  ( 


Svq\-2.5, 
Ks  ' ' 


(17A) 


where  Ks  is  the  equivalent  sand  grain  diameter  and  is  used  as  a measure 
of  the  foil  roughness.  The  larger  value  of  as  predicted  from  equations 
16A  and  17A  is  used. 
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Values  of  Ks  from  References  4 and  5 for  different  surfaces  are  presented 
in  the  following  table. 


SURFACE 


EQUIVALENT  SAND  GRAIN 
DIAMETER  Ks  INCHES[2.54  cm] 


newly  launched  ships 
camouflage  paint  on  airplanes 


.012 

.003 


The  value  of  Ks  that  is  selected  depends  upon  how  the  hydrofoil  is 
used.  Hydrofoils  which  retract  out  of  the  water  and  are  unpainted  will 
have  Ks  values  near  those  given  for  natural  sheet  metal.  Daily 
cleaning  could  even  reduce  these  values.  Hydrofoils  which  are  generally 
kept  submerged  will  probably  have  a Ks  value  in  the  range  given  for  mass 
production  paint  spayed  surfaces.  Fouling  of  the  foil  surface  will 
increase  the  Ks  values  substantially  above  that  of  the  mass  production  spray 
painted  surface. 

D^  represents  the  eddy  making  drag  and  is  given  by  Equation  18A  (See 
Reference  6. 


H 


mi rror 

0 

.00008 

polished  metal  or  wood 

.00008  - 

.00016 

■i 

natural  sheet  metal 

.00016  - 

.00032 

optimum  paint  sprayed 

.00032  - 

.00048 

j 

average  paint  sprayed 

.00048  - 

.00159 

mass  production  paint  sprayed 

.00159  - 

.016 

■5 

i r 


D,  = 0,(2^  60(^)'>) 

avg  avg 
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18A 


where  and  are  the  average  thickness  and  chord  of  the  foil.  The 
above  euqations  for  Dp  D2,  D^,  D^,  Alp  and  AL2  are  all  expressed  in 
terms  of  r^,  the  circulation  strength  at  y=0.  For  the  assumed  elliptic 
circulation,  is  expressed  by  equations  19A,  20A,  and  21A. 


'0  = 


r 


U*c.*a,*a„ 
0 a e 


1+1  0 . 


f2(X,0)  + - Og  f-J  (X,0) 


[93a]  (19A) 


where 


f^(x,a)= 


and 


b'  — 1 2B 

'Vl+4X^+2X 

IVi+4a^/  , 

1 

4xVl+4x^ 

f 

'V^(1+4X^)^/'^ 

1 (x  + log—) 

“na  a 


"0(?) 


^2^^  ,o)“ 


Vl+4X^ 


1 + 0(1 , 1 ) 


[93b]  (20A) 


[93c]  (21A) 


The  above  equations  for  foil  lift  and  foil  drag  are  used  in  the  FOIL 
subroutine  of  the  TRIBRID  computer  program.  They  are  subject  to  the  following 
limitations: 

a)  Fully  wetted,  non  cavitating  flow  exists  over  the  foil  surface. 

This  limitation  is  inherent  in  the  linearized  lifting  line  theory. 

b)  The  foil  has  elliptical  planform.  This  limitation  is  the  result 
of  assuming  an  elliptical  distribution  of  the  circulation  strength. 
Reference  3 states  that  for  an  elliptical  distribution  of  r(y), 
the  planform  is  also  nearly  elliptical  and  that  the  sum  of  the  foil 
drag  components  Dp  D2  + D^,  and  D^  is  approximately  minimum 
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for  a given  lift.  Reference  4 also  recognizes  that  for  a fixed 
geometrical  form  the  elliptic  circ.lation  distribution  can  not 
be  held  constant  for  different  values  of  h.  However,  the 
calculation  of  the  induced  downwash  velocity  at  the  hydrofoil 
supports  the  assumption  that  the  same  circulation  distribution 
will  hold  for  a wide  range  of  depths. 

c)  Foil  aspect  ratio,  AR,  is  not  to  be  lower  that  4.0.  It  is 
assumed  that  the  aspect  ratio  of  the  hydrofoil  is  so  large 
that  the  whole  hydrofoil  may  be  replaced  by  a lifting  line. 
However,  experience  with  airfoils  indicates  that  the  lifting  line 
theory  can  be  expected  to  be  sufficiently  accurate  for  foils  of 
aspect  ratio  AR  as  low  as  4.0. 

d)  The  foil  prediction  method  is  valid  for  only  a certain  range  of 
operating  depths.  The  depth-half  span  ratio  \ should  be  less 
than  1.0.  Equation  6A  is  an  approximation  for  aL^  which  was 
derived  for  x->0.  For  values  of  X greater  than  1 a different 
equation  for  aL^  results,  which  is  not  presented  here  because 

in  the  TRIBRIO  program  we  are  not  interested  in  foil  performance 
at  depths  greater  than  the  half  span  of  the  foil.  The  requirement 
fo>'  fully  wetted  foil  invalidates  the  performance  prediction  at 
zero  depth  or  very  small  depths  where  the  foil  might  be  planing. 
Since  no  prediciton  technique  for  the  minimum  foil  depth  at  which 
the  foil  remains  wetted  is  presented,  the  establishment  of  this 
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minimum  is  a matter  of  engineering  judgement  and  experience. 
Calculations  for  foil  depth  to  average  chord  ratios  as  low  as 
1/2  yield  reasonable  results  and  it  is  suggested  that  this 
value  be  used. for  minimum  foil  depth. 


SYMBOLS  FOR  APPENDIX  B 


A, . . .A 
AF 
ARC 

B,  . . . B 

BSH 

CD 

^F 

Cl 

CP 

CWP 

D 

areo 


n 


n 


D. 


^HE 


^HF 


D 


HSS 


RA 

DSF 

DSE 


^tip 


Chebychev  coefficients  for  hull  and  strut 
frontal  area  of  ship 

arc  length  along  midship  section  of  strut  (see  Fig  IB) 

Chebychev  coefficients  for  hull  and  strut 

maximum  thickness  at  top  of  strut 

aerodynamic  drag  coefficient 

friction  coefficient 

I 3 

strut  waterplane  inertia  coefficient  = 

hull  prismatic  coefficient 

strut  waterplane  area  coefficient 

aerodynamic  drag  of  the  struts  and  superstructure 
above  water 

horizontal  hull  diameter 
hull  eddy  making  drag 
hull  friction  drag 

total  drag  of  the  hulls,  struts  and  superstructure 

drag  due  to  roughness  allowance 

strut  friction  drag 

strut  eddy  making  drag 

tip  drag  for  strut  with  no  hull  underneath 


► 
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DV 

'w 

EXP 

HS 

HWL 


HWS 

'yy 

LCBOLH 


LCFOLS 

LH 

LS 

^SB 

SWS 

^max 

TZ 

U 

XWL 

g 

X. 


vertical  hull  diameter 
wave  making  drag 

strut  flare  exponent 
strut  height  (see  fig  IB) 

strut  depth  below  static  water  line  (ship  is  100% 
buoyant  supported) 

hull  wetted  surface 

strut  waterplane  longitudinal  inertia  about  midship 

ratio  of  hull  center  of  buoyancy  measured  from  midship 
to  hull  length 

ratio  of  strut  center  of  flotation  measured  from 
midship  to  strut  length 

length  of  hull 

hull  buoyant  lift 

length  of  strut 
strut  buoyant  lift 

strut  wetted  surface 

maximum  thickness  of  strut  at  the  waterline 

midship  section  strut  thickness  at  strut  bottom 
free  stream  velocity 
strut  immersed  depth  at  speed 
gravitational  constant 

longitudinal  strut  position  forward  of  midship  section 
longitudinal  hull  position  forward  of  midship  section 
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•strut  half  breadth  at  any  fore  and  aft  location 


hull  half  breadth  at  any  fore  and  aft  location 
measured  in  the  waterplane  through  the  hull  axis 

density  of  air 

density  of  water 

Correlation  Coefficient 
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The  buoyant  lift  provided  by  each  hull  and  each  strut  is  given  in 
Equation  IB  and  2B.  This  buoyant  lift  is  based  upon  the  displacement 
below  an  idealized  level  waterline. 


L^g  = (7r/4)(DV)(DH)(LH)(CP)pg  (IB) 

4b  = pg[(TZ)(XWL)+(||^)(^i^^^^  (2B) 


The  equation  for  43  is  modified  for  air  cushion  supported  concepts  in 
order  to  account  for  the  depressed  waterline  on  the  cushion  side  of  the 
strut. 


FIG.  IB  SINGLE  HULL-STRUT-SUPERSTRUCTURE  GEOMETRY' 
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The  drag  of  the  hull -strut-superstructure  combination  is  given  by  Equation 
3B 


(3B) 


The  individual  component  drags  are  the  following,  Hull  friction  drag 
for  each  hull  is: 


D^P  = Cp-l/2-pU^-HWS  (4B) 

The  hull  wetted  surface  HWS  is  determined  by  dividing  the  hull  into  20 
segments,  calculating  the  average  circumference  and  slant  height  of  each 
segment  and  summing  the  wetted  surface  of  each  segment.  The  resulting 
wetted  surface  area  is  then  corrected  by  subtracting  the  area  masked  by 
the  strut. 

The  shape  of  the  hull  is  generated  by  a series  of  Chebychev  coefficients. 
The  equations  for  obtaining  the  hull  offsets  are  presented  under  the 
discussion  of  wave  drag  since  it  is  the  wave  drag  calculation  which 
requires  that  the  hull  be  represented  by  Chebychev  coefficients. 

The  strut  friction  drag  is  given  by 

D^P  = Cp-1/2P-  U^-SWS  (5B) 

and  the  strut  wetted  surface  SWS  is  given  by 

SWS  = 2((XWL)(LS)+1.33(ARC  - XWL))  (6B) 
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For  cushion  supported  concepts,  SWS  is  corrected  for  the  depressed  waterline 
inside  the  cushion  and  for  struts  with  wetted  bottoms  an  additional  area 
equal  to  (CWP)(LS)(TZ)  is  included.  The  hull  eddy  making  drag  is  given 
by 

= 0.0001  p/2  HWS  (7B) 

and  it  is  based  on  the  residuary  resistance  of  deeply  submerged  bodies  of 
revolution  with  length  to  diameter  ratios  of  10.0  or  greater. 

The  strut  eddy  making  drag  (Reference  3)  is  given  by 

OsE  = + 100(^)'^)SFD 


tip 


added  to  the  eddy  making  drag. 

Tz3 

'tip  TT 


Ty3 

D..._  = 0.0315  ^ CWP  p 


(9B) 


1 2 

(8B) 

■i 

I i 

is 

1 

■J 


The  roughness  allowance  is  given  by 

The  sum  of  the  wave  .making  drag  components,  zD^^,  is  calculated  by  a 
computer  program  written  by  Lin  (Reference  7).  This  program  is  based  on 
linearized  thin  ship  theory  and  assumes  incompressible,  inviscid  potential  flow. 
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The  theory  also  requires  that  the  ship  either  have  a small  camber  which 
satisfies  the  no-cross-flow  condition  between  the  SWATH  type  hulls  and 
struts  or  that  the  thickness  of  each  strut  be  small  compared  to  the 
separation  between  the  struts. 

The  wavemaking  interference  drag  components  are  due  to  the  inter- 
ference of  wave  trains  formed  by  the  different  hulls  and  struts. 

Figure  2B  shows  the  wavemaking  drag  components  that  are  calculated 
for  three  different  basic  configurations.  It  should  be  noted  that 
the  tripple  strut-single  hull  configuration  should  have  an  additional 
wavemaking  drag  component  included  to  account  for  the  wavemaking 
interference  between  the  two  outboard  struts.  This  component  is  expected 
to  be  siiiall  due  to  the  large  distance  between  the  struts  and  to  their 
relatively  small  displaced  volume  at  speed  when  they  are  expected  to  be 
mostly  out  of  the  water. 


! Wave  drag  of  one  strut 
Wave  drag  of  one  hull 
3.  Hull-strut  interference  wave  drag 
•1.  Strut  cross  interference  wave  drag 
r>.  Hull  cross  interference  wave  drag 
6.  Hull -strut  cross  interference  wave  drag 


FIG.  28  SCHEMATIC  DIAGRAM  OF 

WAVE  DRAG  AND  WAVE  INTERFERENCE  DRAG  BETWEEN  HULLS  Af;0  STRUTS 
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The  computer  program  for  calculating  the  wave  drag  requires  that 
the  hull  and  strut  shape  be  defined  by  a series  of  Chebychev  coefficients. 
With  many  Chebychev  coefficients  (up  to  15  have  been  used)  almost  any 


desired  hull  and  strut  shape  can  be  generated  to  within  1%  accuracy, 
however,  in  order  to  save  computer  time,  only  3 Chebychev  coefficients 
are  used  in  this  report. 

The  hull  Chebychev  coefficients  A and  B are: 


= 4-CP/PI 
A2=1-A, 

= 2-A^-LCBOLH 
^3  B2=A2=0.0 


(IIB) 

(12B) 

(13B) 

(14B) 


and  the  strut  Chebychev  coefficients  are: 


4-CWP/PI 

(15B) 

Aj  = 

^ - 

64-CI/PI 

(16B) 

*3  ' 

1 - 

Ai  - A^ 

(17B) 

2-A, 

•LCFOLS 

(18B) 

=2' 

= 

0 

(19B) 

The  equation  for  calculating  the  offsets  of  the  strut  at  any  waterline  is 

yg(x)g  = + y^gf  [A^.cos[(2i-l)e]  + B.  sin[2ie]]  (20B) 

■^ref  breadth  of  the  waterline  at  midships  and  is  constant 

for  struts  with  no  flare. 


52 


The  offset  of  the  hull  at  a waterline  through  the  axis  of  rotation 
of  the  hull  is 

cos[(2i-l)e]  + s1n[2le]]  (22B) 


-1 

where  8=  sin  {-^) 


(23B) 


Xs 


FIG.  3B  STRUT  OFFSETS  (plan  view  of  strut-waterlihe) 

I 


DH 

For  hull  offsets,  replace  subscript  sbyh  and  y^^^  by  -y  in  Figure  3B. 

The  aerodynamic  drag  component  for  the  parts  of  the  ship  above 

water  is  given  by 


0 


areo 


CD-p  -U^-AF/2 

a 


(24B) 


Reference  6,  p.  11-8  presents  values  of  CD  = .68  for  a passenger  liner, 
1.22  for  a cargo  ship  and  0.17  for  a hypothetical  streamlined  ship. 

A CD  value  between  0.17  and  1.22  is  suggested  depending  on  the  extent  of 
streamlining. 


w- 
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LIST  OF  SYMBOLS  FOR  APPENDIX  C 


cuss 


seal 


TS 


wp 


D 


mom 


‘"seal 

FB 

FIP 

F. 


GRAV.g 
LC,  1 
LS 

^stat 

Pc 

R 

SEAL  W 
SEP 

^CWL 


cushion  area 

total  wetted  area  of  the  seals 
friction  coefficient  for  seal 
waterplane  area  coefficient  of  each  strut 
momentum  drag 
seal  drag 

Froude  number  based  on  cushion 

Fan  ideal  power 

Froude  number  based  on  £ 

Froude  number  corresponding  to  hump  speed 

gravitational  constant 
cushion  length 
strut  length 
aerostatic  lift 

cushion  pressure 

wave  resistance  of  pressure  patch 
width  of  seal 

separation  distance  between  strut  centerlines 

maximum  strut  thickness  at  the  elevation  corresponding 
to  the  waterline  on  the  cushion  side  of  the  strut. 
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u 

WPAR 

b 


velocity  of  pressure  patch 

wave  resistance  parameter 

beam  of  pressure  distribution 

Eulers  constant  .577... 

length  to  beam  ratio  of  pressure  patch 

water  density 

air  density 


AEROSTATIC  LIFT 


The  lift  provided  by  the  air  cushion  is  given  by: 
^STAT  " ^ \uss 


(1C) 


the  cushion  area  is  estimated  by  equation  2C  which  is  exact  for  any  value  of  LC 
provided  that  CUP  = Ij and  for  any  value  of  CWP  provided  that  LC  = LS 

\uss  = * SEP  - C„p  ♦ LC  * (20 

The  air  flow  required  to  maintain  cushion  pressure  is  given  by  equations  3C 

and  3D.  Equation  3C  is  based  on  model  experimental  results  and  was  developed 
by  Mr.  Robert  A.  Wilson  of  DTNSRDC  (Code  1630). 

FB  = U/  / GRAV  * SEALW  (3C) 

Q = 1.67  (0.203  FB^  + .15  FB  + .247)  x SEALW^  * /F  (4C) 

The  ideal  lift  power  is  FIP  = Q * P^ 
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DRAG  DUE  TO  AEROSTATIC  LIFT 


The  wavemaking  resistance  due  to  the  moving  cushion  pressure  patch 
is  estimated  using  the  dosed  form  approximations  derived  by  Ravenscroft 
in  Reference  8 . The  wave  resistance  parameter  WPAR  is  given  by  equation 


5C,  6C  or  7C  depending  on  the  Froude  number  F^. 


4 * i/4  2 ^ 3 - Y + 

WPAR  = ,rF„  (1+1.6 


d + 1.6  F„  ^ 


16X 


(5C) 


-/X  F o 
WPAR  = .2 

e sin 


u. 


1 


{ 3 - Y 


+ J.n 


^ ^ , 1 


4 ,1/4 


4X 


} <F  <fx--F 


(6C) 


WPAR  = e 


-/X  F r 

^ . 2 
sin 


2 F. 


■+  i 

4tt 


^ {3  - Y + Jin  (-|-)}  ; 


4 

I’- 


(7C) 
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4 


The  primary  hump  F.Poude  number  F,  is  estimated  by  Equation  8C 

•H 


F„  =\A+  .01745x2  . 

V, 


(8C) 


The  wave  resistance  due  to  the  moving  pressure  patch  is  given 


by 


R = WPAR  X 4 P b/(  *GRAV). 


(9C) 


The  momentum  drag  associated  with  the  lift  system  is  simply  given  by 


D = U*p  . *0 
mom  "^air  ^ 


(IOC) 


where  Q is  given  by  Equation  4C  and  the  seal  drag  based  on  flat  plate 
friction  is  given  by 


0 • “ C__ • p • U • A 1 

seal  FS  seal 


(nc) 
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(1)  DTNSRDC  REPORTS.  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE.  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

!2)  OEP.ARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES.  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE.  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
REPORTS;  NOT  FOR  GENERAL  DISTRIBUTION. 


